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ABSTRACT 

Let  X^,  . . . ,  X^  and  . be  two  independent 

random  samples  from  populations  with  continuous  IFRA  distribu¬ 
tions  F(-)  and  F(*/A)  respectively,  and  let  s^, 
denote  the  ranks  of  the  Y’s  in  the  combined  sample.  For 
testing  Hq:  A  £  1  vs.  H^:  A  >  1  it  is  shown  that  the 
error  probabilities  of  each  monotone  rank  test  $  are  bounded 
by  the  error  probabilities  for  exponential  alternatives,  i.e. 
if  K(t)  -  1  -  exp(-t),  t  ^  0,  then  P(reject  Hq|f(*)»  F(-/A)) 
<:  P (reject  Hq|k(*),  K(./A))  for  A  _<  1,  and  P(accept 
Hq | F C * ) .  F( • /A) )  £  P(accept  H  |K(*)»  K(-/A))  for  A  >  1  . 

These  inequalities  are  used  to  derive  tests  that  maximize  the 
minimum  power  in  the  class  of  all  rank  tests.  The  results 
are  extended  to  censored  samples,  sequential  sampling,  dis¬ 
tributions  ordered  by  skewness,  the  problem  of  combining 
Independent  test  statistics,  and  the  goodness  of  fit  problem. 


MINIMAX  RESULTS  FOR  IFRA  SCALE  ALTERNATIVES 

by 

w 

Kjell  Doksum 

1.  INTRODUCTION  AND  SUMMARY 

» 

Recent  results  by  Birnbaum,  Esary  and  Marshall  (1966),  Barlow  and  Proschan 
(1967)  and  others  suggest  that  the  exponential  models  used  by  Epstein  and  Sobel 
(1953)  and  others  for  life  testing  problems  should  be  extended  to  models  in  which 
the  lifetimes  have  increasing  failure  rate  average  (IFRA)  distributions. 

In  this  paper,  IFRA  scale  models  are  dealt  with.  For  instance,  consider  two 
independent  random  samples  X^,  ...»  and  Y^,  ...»  Y^  from  populations  with 

distributions  F(*)  and  F(-/A)  and  respectively,  where  F(*)  is  a  continuous, 
unknown,  IFRA  distribution.  The  null  hypothesis  H^:  A  <_  1  is  to  be  tested  against 
the  alternative  H^:  A  >  1  .  Let  K(-)  be  the  exponential  distribution  function 

defined  by  K(t)  *=  1  -  exp(-t),  t  _>  0  .  It  is  shown  (Theorem  2.1)  that  for  each 

monotone  rank  test  and  each  A  >  1,  the  minimum  power  is  attained  at  the  ex¬ 
ponential  distribution,  i.e.  inf^#j  E[$|F(‘),  F(’/A)]  *  E[$|K(*)»  K(  *  /A)  ]  . 

Moreover,  it  is  shown  that  for  each  A  <  1,  the  probability  of  a  type  I  error  is 
maximized  at  the  exponential  distribution,  i.e.  sup^,^  E[$|F(*)»  F( •  /A )  ]  * 

-  E[*|K(*),  K(-/A)]  . 

Since  F()  is  unknown,  it  is  not  possible  to  maximize  the  power 

E[$|F(')»  F ( * /A)  ] ,  A  >  1  .  However,  it  is  shown  (Theorems  3.1,  3.2,  4.1,  4.2)  that 

the  tests  that  maximize  the  power  for  the  exponential  alternative  K(’/A)  actually 
maximize  the  minimum  power  infj,^  E[$|f(*),  F  ( *  /A )  ]  .  Thus  these  tests  are  minimax. 
They  have  been  computed  by  Lehmann  (1953),  Savage  (1956),  and  Rao,  Savage  and  Sobel 
(19^0).  The  results  indicate  that  in  the  case  of  uncensored  samples,  one  should 
use  one  of  the  statistics 
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(1.1) 

(1.2) 


L  =  l  log  [1  -  s  /(N+j)],  or 

J-l  J 


n 

S  =  V  J  (s  ) ,  with  J  (k) 

j-l  °  j  0 


n 


I 

j=N+l-k 


1/J  . 


where  N  -  m  +  n  and  s,  ,  ....  s  are  the  ordered  ranks  of  the  Y's  in  the 

1  n 

combined  sample  of  X's  and  Y’s  .  L  is  minimax  for  A  in  an  interval  about 
two,  and  S  is  minimax  for  A  in  an  interval  (1,A)  to  the  right  of  one. 

The  minimax  statistics  in  the  case  of  censored  samples  are  more  complicated 
(see  (4.1)  and  (4.2))  and  one  might  use  one  of  the  approximations  suggested  by 
Gastwirth  (1965)  or  Basu  (1967)  (see  (4.3)). 

Next,  bounds  on  the  error  probabilities  for  sequential  rank  tests  are  given 
in  Section  5,  and  inequalities  showing  that  power  is  a  decreasing  function  of 
"skewness  to  the  right"  are  given  in  Section  6.  Finally,  extensions  to  the  problem 
of  combining  independent  tests  statistics  are  given  in  Section  7. 

Only  finite  sample  size  properties  are  dealt  with.  Asymptotic  results  are 


given  in  [8] . 
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2j. _ TH1.  POWER  OF  MONOTONE  TESTS 

Let  X^,  . ..,  and  Y^ . Y^  be  two  independent  random  samples  from 

populations  with  distribution  functions  F(*)  and  G(*)  .  The  problem  first  con¬ 
sidered  is  that  of  testing  the  null  hypothesis  that  the  X's  are  stochastically 
larger  than  or  equal  to  the  Y's  against  the  alternative  that  the  X's  are 

stochastically  smaller  than  the  Y's;  i.e.  one  tests  H  :"G(t)  >  F(t)  for  all  t," 

o 

against  H^:"G(t)  <_  F(t)  for  all  t  and  G(t)  <  F(t)  for  some  t,"  <  t  <  00  . 

For  this  problem,  it  is  reasonable  to  consider  tests  $  that  are  monotone  in 
the  sense  that 

■ 

<_  y^  for  j  ■  1,  ...»  n  implies 

(2.1)  *(xl . V  yl . V 

i  *(*i . v  yi . V 


Monotone  tests  are  known  (e.g.  [12,  p.  187]  and  [4])  to  have  monotone  power,  i.e. 
if  6 («t> ; F ,G)  =  E [ 4> |  (F , G)  ]  denotes  the  power  of  then  G^(‘)  <_  G2(*)  implies 

8 (<J> ; G2  » E)  <_  8(^;C^,F)  .  In  particular, 

Lemma  2 . 1 

If  4>  is  monotone,  then 


(2.2) 


(2.3) 


G^(t)  <_  G2(t)  _<  F(t)  for  all  t  implies 
8(4> ;F,F)  <_  8(<f>;F,G2)  <_  6 (<J> ; F ,G1>  and 
F(t)  <_  G^(t)  <_  G2(t)  for  all  t  implies 
8($;F,G2)  £  B(<f;F,G1)  £  B($;F,F) 
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This  result  shews  that  for  monotone  tests,  the  further  the  Y's  are  away  from  the 
X's  stochastically,  the  smaller  are  the  probabilities  of  type  I  and  type  II 
errors. 

Next  the  scale  model  in  which  G(t)  =  F(t/A)  for  all  t  and  some  A  >  0 
will  be  considered.  The  problem  is  now  to  test  "A  <_  1"  against  "A  >  1." 

For  each  test  <J>,  let  8(<J>;F;A)  «  E[$|f(*)»  F  ( •  /  A) )  denote  the  power  of  $  for 
the  scale  alternative.  F  •=  F  ( * )  is  said  to  be  an  IFRA  distribution  if  F(0)  *  0 
and 

(2.4)  -log[l  -  F ( t )  ] /  is  r.ondecreasing  in  t  >  0  . 

Similarly,  F  is  a  DFRA  distribution  if  F (0)  *=  0  and  -log[l  -  F(t)]/t  is  non¬ 
increasing  in  t  >  0  .  Let  s,  <  s„  <  ...  <  s  denote  the  ordered  ranks  of 

12  n 

Y^ . Y^  in  the  combined  sample  and  let  rank  tests  be  tests  depending  only  on 

81»  *  *  *  »  sn  * 

It  will  now  be  shown  that  for  monotone  rank  tests,  the  probabilities  of  type 
I  and  type  II  errors  for  IFRA  distributions  are  at  most  what  they  are  for  the 
exponential  distribution,  i.e.  if  K(t)  =  1  -  exp(-t),  t  _>  0,  then 

Theorem  2.1 

If  4>  is  a  monotone  rank  test  and  F  is  a  continuous  IFRA  distribution, 

then 

for  each  A  _>  1,  3(4>;K;A)  <_  B ( ; F ; A)  ,  and 
for  each  0  <  A  <_  1,  B  ( <p ;  F ;  A)  <_  BUjKjA)  . 


(2.5) 

(2.6) 

Proof : 

Since 


B(t)  =  -log[l  -  F ( t ) ]  is  an  increasing  function  and  since  the  ranks 
are  invariant  under  increasing  transformations,  then 


5 


(2.7) 


*(X1 . V  Y1 . Yn) 

-  ♦(X1,  ....  Xm;  Yl,  Yn) 


where  *=  B(Xi>  and  «=  B(Y^)  . 


Next  note  that  if  one  defines  Y^  -  -A  log|l 
IFRA  property  implies  that  for  each  A  _>  1, 


\  A/J’ 


then  the  definition  of  the 


(2.8) 


i  i  < 

log[l  -  F(Y  /&)]  <  -  ~  log[l  -  F(Y  )]  -  y1  . 

J  J  J 


f  f  t 


That  is,  Yj  £  Y ^ ,  and  sinv_e  $  is  monotone,  then 


(2.9) 


*(X1 . Xm’  Y1  '  •**»  Yn  J 


S.  ^(^i»  •••»  X  ;  Y.  ,  . .  •  |  Y  ) 
x  ml  n 


Upon  using  (2.7),  taking  expectations  in  (2.9),  and  noting  that  has  the  dis- 

1 1 

tribution  K(*)  while  Y^  has  the  distribution  K('/A),  one  obtains  (2.5).  The 
inequalities  are  reversed  for  0  <  A  £  1,  and  (2.6)  follows. 

For  DFRA  distributions,  the  opposite  result  holds,  i.e. 

Lemma  2.2 

If  4>  is  a  monotone  rank  test  and  F  is  a  continuous  DFRA  distribution,  then 


(2.10) 


for  each  A  >_  1,  B($;F;A)  <_  0($;K;A),  and 


(2.11) 


for  each  0  <  A  <_  1,  B($;K;A)  <_  8(<f>;F;A) 


The  proof  is  the  same  as  that  of  Theorem  2.1  with  the  inequalities  reversed. 


■  -4a 
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Remark 

Lemma  2.1  shows  that  if  $  is  monotone,  then  B(4>;F;A)  is  an  increasing 
function  of  A  >  0  .  Thus  Theorem  2.1  implies  that  the*  power  of  monotone  rank 
tests  has  the  bound  f}($;F;A)  j>  6(4>;K;A^)  for  the  IFRA  scale  alternative  with 
A  _>  Ax  _>  1  • 
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3.  MINIMAX  TESTS  IN  THE  TWO-SAMPLE  CASE 

Although  It  is  a  desirable  property  for  a  test  to  have  monotone  power  for  the 
stochastic  ordering  alternative  of  Section  2,  this  alternative  cannot  be  used 

to  derive  optimal  tests.  One  would  like  to  have  tests  that  are  most  likely  to  re¬ 
ject  for  alternatives  that  indicate  a  definite  distinction  between  the  distributions 
of  X  and  Y.  When  X  and  Y  are  time  measurements,  then  one  such  model  is  the 

scale  model  of  Section  2  in  which  Y  has  the  same  distribution  as  AX  for  some 

A  >  0  .  Then  [ E (Y )  -  E(X)]/E(X)  -  A  -  1,  and  (A-l)  measures  the  relative  differ¬ 
ence  of  the  mean  times. 

I 

Note  that  if  N  *=  m  +  n  and  if  s^  is  detined  to  be  the  number  of  Y's  greater 

than  or  equal  to  the  (N+l-j)th  order  statistic  in  the  combined  sample,  then 
»  * 

(s^,  ...,  sN)  is  equivalent  to  the  ordered  ranks  (s^,  ...,  s^)  .  Next  it  will  be 

shown  that  in  the  class  of  level  a  rank  tests,  the  Lehmann  level  a  test 

which  rejects  for  small  values  of 


maximizes  the  minimum  power  inf  0 (<h  * F; A)  for  the  IFRA  scale  model  whenever  A  >  1  . 
Here,  a  rank  test  <f>  is  said  to  be  of  level  a  if  6 (4> ; F ; F)  *  E(<j>|(F,F))  *  a  for 
all  continuous  distribution  functions  F  .  As  in  Section  2,  0(4>;F;A)  * 

*  E[4»  |F(  * )  ,  F(  *  /A)  )  denotes  the  power  of  <f>  for  the  scale  alternative. 

Theorem  3.1 

The  Lehmann  test  ^  is  minimax  for  the  IFRA  scale  alternative  in  the  sense 
that  for  each  scale  parameter  A  >  1  and  for  F  ranging  over  the  class  of  con¬ 
tinuous  IFRA  distributions } 

a 

(3.2)  inf B($;F;A)  <  infp  B(^;F;A) 


- 
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for  each  level  a  rank  test  4>  .  Moreover t  ^  ia  the  unique  minimax  rank  test  • 

in  the  sense  that  any  other  level  a  rank  test  satisfying  (3.2)  coincides  with 

a.e. 

* 

Proof : 

— 

Using  the  Neyman-Pearson  Lemma,  Lehmann  (1953),  Savage  (1956),  and  Rao,  Savage 
and  Sobel  (1960,  Corollary  3.4)  have  essentially  shown  that  $ maximizes  the 
power  for  the  exponential  scale  model,  i.e.  for  each  A  >  1, 

(3.3)  P($;K;A)  £  8(^;K;A) 

i 

r 

for  all  level  a  rank  tests  $  .  On  the  other  hand,  the  definition  (3.1)  shows 
that  <f>A  is  monotone,  thus  Theorem  (2.1)  implies  that  attains  its  minimum 

power  for  the  exponential  scale  model,  i.e. 

♦ 

(3-4)  infp  B($a;F;A)  -  B(^;K;A) 

Since  K  is  an  IFRA  distribution,  then 

(3.5)  infp  B(<j>;F;A)  <  B($;K;A) 

and  (3.2)  follows.  Uniqueness  follows  also  since  if  ip  is  a  level  a  rank  test 
satisfying  (3.2),  then  the  above  shows  that  it  must  satisfy 

(3.6)  B(*;K;A)  £  B(*;K;A) 

for  all  level  a  rank  tests  <J>  .  Since  ^  also  satisfies  this  (see  (3.3)),  then 

the  uniqueness  part  of  the  Neyman-Pearson  Lemma  implies  that  <Ji  *  a.e.  ^ 

In  order  to  be  able  to  use  the  minimax  test  one  must  choose  a  value  of  A  . 

Savage  (1956)  suggests  using  the  level  a  test  $  that  maximizes  the  power  for 

8 

the  exponential  model  for  A  in  a  neighborhood  (1,6)  to  the  right  of  one,  i.e. 


9 


when  the  relative  difference  of  the  means  of  X  and  Y  is  close  to  zero  (and 
positive).  This  test  reject  for  large  values  of 


(3.7) 


l  J  (s  ) ,  where  J  (k)  -  ][  1/j  . 

j-1  J  J-N+l-k 


Note  that  Jo(k)  can  be  approximated  by  -log[l-k/(N+l)  ]  .  It  can  now  be  shown 

that  the  Savage  test  $  is  minimax  for  A  in  an  interval  of  the  form  (1,6)  . 

s 

Theorem  3.2 

The  Savage  teat  $  ie  locally  minimax  for  the  IFRA  scale  alternative  in  the 
sense  that  there  exists  6  >  1  such  that  for  F  ranging  over  the  class  of  con¬ 
tinuous  IFRA  distributions , 


(3.8) 


inf p  0(«P;F;A)  £  infp  B(*8;F;A) 


for  each  level  a  rank  test  $  and  for  all  A  in  the  interval  (1,6)  .  More¬ 
over,  <p  is  the  unique  locally  minimax  rank  test  in  the  sense  that  any  other 

8 

level  a  rank  test  satisfying  (3.8)  coincides  with  $  a.e. 


Proof : 


Savage  (1956)  has  essentially  shown  that  there  exists  6  >  1  such  that 


(3.9) 


0(4>;K;A)  £  B(*  ;K;A) 


for  all  level  a  tests  $  and  all  A  in  (1,6)  .  Since  p  is  monotone,  the 
remainder  of  the  proof  is  as  the  proof  of  Theorem  1. 


Remarks 


(i)  The  Savage  test  maximizes  the  minimum  power  when  the  relative  difference 
of  the  means  of  X  and  Y  (A-l)  is  close  to  zero.  In  most  situations,  it  would 
be  better  to  use  the  test  that  is  optimal  when  (A-l)  is  in  a  neighborhood  of  some 


T 


1 
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fixed  positive  number  X  (say).  Thus  one  would  use  the  Lehmann  test  ^  defined 
by  (3.1)  with  A  =  X  +  1  .  For  instance,  if  the  relative  difference  of  the  means 
of  X  and  Y  is  unity,  then  Lehmann  (1953)  has  shown  that  ■  $2  i®  equiva¬ 
lent  to  the  test  that  rejects  for  small  values  of 

(3.10) 


jli <H  +  J  -  V  or  j,  M1-  h) 


(ii)  The  uniqueness  parts  of  Theorem  3.1  and  3.2  can  be  extended  as  follows. 

A  test  <f>  is  said  to  be  distribution- free  (DF)  if  B(<j>;F,F)  -  E($|(F,F))  is 

independent  of  F  for  F  continuous.  Thus  rank  tests  are  DF.  For  A  close  to 

one,  tests  that  are  not  DF  have  minimum  power  less  than  $  and  $  ,  and  they 

cannot  be  mlnimax.  To  see  this,  note  that  if  41  is  of  level  a,  a.e.  continuous, 

and  not  DF,  then  there  exists  a  continuous  distribution  F  such  that 

o 

E[$|f  ,F  )]  <  a  -  sup  E [ d)  |  (F ,F)  ]  .  Now  for  A  close  to  one  B($;F  ;A)<  a  and  <p 
o  o  r  o 

is  worse  than  d>,  and  A 
A  s 

(iii)  Lemma  2.1  shows  that  if  $  is  a  monotone  test,  then  8($;F;A)  is  an 

increasing  function  of  the  scale  parameter  A  >  0  .  This  implies  that  for  fixed 

Aj>  1,  $  is  doubly  miniraax  for  the  scale  alternative  in  the  sense  that  for 

testing  H  :  A  £  1  against  H  :  A  ^  A.  ,  it  maximizes  inf  [inf.  .  B($;F,A)]  ■ 

^  O  XX  r 

-  inf  .  [inf  6 (d> ;F ; A)  ]  .  Here,  F  ranges  over  the  class  of  continuous  IFRA 
distributions  and  only  level  a  rank  tests  are  considered.  See  the  remark  at  the 
end  of  Section  2.  Note  that  (1,A^)  is  an  indifference  region. 
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4.  MINIMAX  TESTS  BASED  ON  CENSORED  SAMPLES 


The  censored  samples  considered  here  arise  typically  as  follows:  m  objects 

of  one  type  and  n  objects  of  a  second  type  are  put  on  trial  at  the  same  time  and 

*  * 

one  waits  until  a  total  of  N  <  m  +  n  -  N  of  the  objects  have  failed,  where  N 

is  a  fixed  number.  Moreover  the  experiment  is  conducted  so  that  if  one  waited 

until  all  N  ■=  m  +  n  objects  failed,  then  the  times  to  failure  X,,  ...,  X  and 

1  m 

Yj,  ...,  would  be  two  independent  random  samples. 

* 

Thus  the  situation  is  as  in  Sections  2  and  3  except  that  only  the  first  N 

A  A 

smallest  order  statistics  in  the  combined  sample  are  observed.  Let  m  and  n 
denote  the  total  number  of  X's  and  Y's  observed  respectively.  Since  the  unob¬ 
served  X's  and  Y's  are  all  larger  than  the  observed  ones,  it  is  possible  to 

compute  the  ordered  ranks  s^  <  ...  <  s^*  of  Y^ . Y^*  among  the  urzcensored 

combined  sample  X^,  ...,  X^;  Y^,  . . . ,Y^  .  Rao,  Savage  and  Sobel  (1960,  Corollary 

A 

3.4)  have  computed  the  most  powerful  test  based  on  s^,  ...,  s^*  for  ex¬ 
ponential  alternatives.  This  level  a  test  rejects  for  small  values  of 


(4.1) 


^A 

l  log["A(A)  +  j  -  8  +  n  log  A 

j-1  L  J 

with  A(A)  ■  (m-m)  +  (n-n  )/A 


where  s^  is  as  in  Section  3.  For  tests  depending  only  on  s^,  ...,  s^*  one  obtains 


Theorem  4.1 


The  test  ^  is  uniquely  minimax  for  the  IFRA  scale  alternative  in  the  sense 
of  Theorem  3.1. 


Proof : 


It  is  enough  to  show  that  is  monotone  and  then  apply  the  proof  of  Theorem 
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However,  It  is  clear  from  the  definition  (4.1)  that  $  is  monotone  since  A(A) 

6 

•  * 

-s  and  n  are  decreasing  functions  of  the  y  s  . 


i 


The  locally  most  powerful  test  $  (Rao,  Savage  and  Sobel  (1966,  Corollary 

8 


3.4))  for  the  exponential  alternative  K(*/&)  rejects  for  large  values  of 


(4.2) 


5  *  v  ,  it  it  it  if 

I  J  (s  )  +  (n-n  )J  (N  )  -  n 

j-1  3 

* 

*  N  *  * 

where  J  (k)  ■  ^  l/[j  +  (m-m  )  +  (n-n  )] 

°  j-N  -k+1 


Theorem  4.2 

The  test  <p  is  uniquely  locally  minimax  for  the  IFRA  scale  alternative  in  the 
s 

sense  of  Theorem  3.2. 

Proof : 

Again,  it  is  enough  to  show  that  4>  Is  monotone.  Monoticity  follows  at  once 

8 

if  it  it  if 

from  (4.2)  and  the  inequality  J  (s.)  <  J  (N  ) ,  j-1,  ...»  n 

o  j  —  o 

Remark 


The  tests  and  are  complicated  and  one  may  use  the  approximations 

suggested  by  Gastwirth  (1965)  and  Basu  (1967).  The  latters  paper  contains  tables 
of  rejection  points  and  power  for  the  test  based  on  the  statistic 

* 

l  Jo(sJ)+!!ZS!f  J*  J  (j)  -  i  (TTtn) 

j-1  2  N-N  j-N  +1 

(A. 3) 


I  J  (s.)  +  (n-n  )  J  (N  )  -  n  +ir  (m-n) 
j-1  °  J  0  2 


N 

where  as  in  Section  3,  J  (k)  -  £  1/j  . 

0  j-N-k+1 
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5.  SEQUENTIAL  TESTS 

Rao,  Savage  and  Sobel  (1960),  Sobel  (1966),  Basu  (1967)  and  others  have  pointed 
out  that  if  one  is  measuring  time  variables,  then  all  rank  tests  have  interpretations 
as  truncated  (curtailed)  sequential  tests.  Consider  for  instance  the  uncensored 
samples  of  Section  3.  If  the  significance  level  a  is  .10,  if  m  *>  n  ■  5  and 
if  the  first  four  values  observed  are  X's,  then  one  can  stop  the  experiment 
since  H  will  be  rejected  no  matter  what  the  remainder  of  the  observations  turns 
out  to  be.  Thus  only  four  observations  are  needed.  The  maximum  number  of  ob¬ 
servations  needed  until  a  decision  is  reached  (the  truncation  point)  is  ten.  Tables 
for  carrying  out  some  truncated  sequential  rank  tests  and  tables  containing  expected 
sample  sizes  and  sample  times  are  given  by  Sobel  (1966)  and  Basu  (1967). 

The  situation  considered  here  is  the  one  in  which  the  observations  at  the  nth 
stage  form  two  independent  random  samples  X.,  ...,  Xr;  Y^,  ...»  Yq  .  A  sequential 
test  $  is  monotone  if 

i 

y^  <_  y^  for  j  -  1,  2,  ...  implies 

(5.i)  4>(xltx2,  y1»y2>  •••) 

t  I 

1  (j»(x1,x2,  ...;  y1.y2*  •••) 


From  the  arguments  of  the  proof  of  Theorem  2.1,  it  follows  that  for  the  IFRA  scale 

model,  the  error  probabilities  of  each  monotone  sequential  rank  test  is  bounded  by 

the  error  probability  for  the  exponential  alternative  K(t/A)  ■  1  -  exp(-t/A)  . 

Here  rank  tests  are  tests  which  at  the  nth  stage  depend  on  the  ordered  ranks 

s, ,  ....  s  of  the  Y's  in  the  combined  sample  X, ,  . . . ,  X  ;  Y,  ,  . . . ,  Y  .  As 

I  n  r  1  n  1  n 

i 

before,  s^  Is  the  number  of  Y's  greater  than  or  equal  to  the  (N+l-j)th  order 

»  * 

statistic  in  the  combined  sample.  Thus  (s^,  ...,  s2n)  is  equivalent  to 

(s.,  ....  s  )  .  0  (<b ;F ; A)  denotes  the  power  of  <f>  when  the  X's  and  Y's  have 

distributions  F( • )  and  F(*/A)  respectively. 


Theorem  5 . 1 


For  each  monotone  sequential  rank  test  $  and  each  continuous  IFRA  distribu- 


tion  F 

(5.2) 

A  <_  1  implies 

0  (d> ;F ; A)  <_  0  U;k;A)  and . 

S  8 

(5.3) 

A  j>  1  implies 

0 .(♦;F;A)  >  0  (4>;K;A) 

o  8 

Proof : 

Same  as  the  proof  of  Theorem  2.1. 

A  A 

The  sequential  probability  ratio  test  $  *  4> (a »b )  based  on  ranks  for  the 
exponential  alternative  K(*/A)  has  been  studied  by  Parent  (1965),  Savage  and 
Sethuraman  (1966)  and  others.  For  testing  A  «*  1  against  A  -  A^  >  1,  it  can  be 
written 


Take  two  more  observations  (X  Y  .,)  if 

n+1 ’  n+1 

log(a)  <_  Ln  <_  log(b)  , 

(5.4)  accept  Hq  if  <  log(a), 

reject  Hq  if  Lr  >  log(b), 
n  ®  1,2 ,  ... 


where  0  <  a  <  1  <  b  are  constants  independent  of  n  and 


(5.5) 


L 

n 


N 

log [ (2n) ! ]  -  n  log  A1  -  [  log  j 

j-l  L 


.6-1) 

A, 


It  is  clear  that  this  test  is  monotone,  thus  the  bounds  on  the  error  prob¬ 
abilities  given  in  Theorem  5.1  apply.  Let  aQ  and  be  desired  error 

probabilities  for  A  ■  1  and  A  ■=  A^  respectively,  and  let  a  and  b  be  such  that 
the  test  <J>(a,b)  achieve  these  bounds  for  the  exponential  distribution,  i.e. 

*  A 

0  (4> ; K ;  1)  ■  a  and  6  (4>;K;A..)  -  a,  .  Then  Theorem  5.1  and  the  sequential  version 
s  o  8  1  1 


15 


of  Lemma  2.1  implies  that  and  a^  are  bounds  on  the  error  probabilities  for 

all  I  FRA  distributions  for  A  £  1  and  A  >_  A^  respectively;  i.e. 


(5.6) 


0  (4>;F;A)  <  a  for  A  <  1,  and 
s  —  o 

A 

8  ( 4> ; F ; A )  >  1-a,  for  A  >  A, 

s  —  l  l 

for  all  continuous  IFRA  distributions  F  . 


Note  that  the  test  $  is  easy  to  carry  out  as  a  and  b  can  be  closely  approx¬ 
imated  as  follows 


(5.7) 


a  1-a 

a  i  -  and  b  -  - - 

1-a  a 

o  o 


Since  $  is  not  known  to  be  optimal  for  the  exponential  alternative  K(*/A), 


it  is  not  possible  to  prove  any  minimax  optimality  properties  of  this  test. 
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6.  GENERALIZATIONS. SKEWNESS 

results  of  the  previous  sections  hold  if  the  class  of  continuous  IFRA 
distributions  is  replaced  by  any  class  of  continuous  distributions  that  is  con¬ 
tained  in  this  class  and  contains  the  exponential  distribution;  e.g.  the  class  of 
continuous  IFR  (increasing  failure  rate)  distributions. 

A  distribution  F(*)  is  said  to  be  IFR  if  K  ^ ( F ( * ) ]  is  convex,  where 
K(t)  *  1  -  exp(-t)  ,  t  _>  0  .  Extending  this,  van  Zwet  (1964)  calls  a  distribution 
F(«)  lees  skew  (to  the  right)  than  H(*)  if  H  ^ [ F ( • )  ]  is  convex,  and  he  writes 
F  <  H  .  Thus  a  distribution  is  IFR  if  it  is  less  skew  than  the  exponential  distri¬ 
bution  K(*)  .  Carrying  these  ideas  further,  we  define  a  distribtuion  F(*)  to  be 
weakly  less  skew  (to  the  right)  than  H  if  H  ^[F(*)l  is  starshaped  (H  *[F(t)]/t 
is  non-decreasing  in  t  >  0) ,  and  write  F  £  H  .  Thus  a  distribution  is  IFRA  if 
it  is  weakly  less  skew  than  the  exponential  distribution  K(-)  .  One  can  now  define 
F  to  be  more  IFRA  than  H  if  F  £  H  .  It  Is  assumed  throughout  this  section 
that  all  the  distributions  satisfy  F(0)  *  0  and  have  inverses  on  the  interval 
(0,°°)  .  From  the  properties  of  starshaped  and  convex  functions  It  follows  that 

Lemma  6 . 1 

(i)  If  F  is  weakly  less  skew  than  H,  and  H  is  weakly  less  skew  than 

L,  then  F  is  weakly  less  skew  than  L  . 

(ii)  If  F  is  less  skew  than  H,  then  F  is  weakly  less  skew  than  H  . 

Next  one  shows  that  for  scale  alternatives,  the  error  probabilities  of  each 
monotone  rank  test  4>  are  decreasing  functions  of  weak  skewness  (to  the  right),  or 
equivalently,  increasing  functions  of  the  IFRA  property. 

Theorem  6 . 1 

If  <p  is  a  monotone  two-sample  rank  test  and  if  F  and  H  are  continuous 
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distributions  such  that  F  is  weakly  less  skew  than  H,  then 

(I)  for  each  A  £  1,  8($;F;A)  £  B (4> ;H ; A) 

(II)  for  each  A  £  1,  8(<$>;H;A)  £  8($;F;A) 

Proof : 

The  proof  is  the  same  as  for  Theorem  2.1  with  H  ^[F(*)l  replacing  -log[l-F(*)] 

One  can  now  proceed  as  in  Section  3  and  show  that  the  test  that  maximizes  the 
power  for  the  alternative  H(*/A)  maximizes  the  minimum  power  inf^  0($;F;A)  over 
the  class  of  F  that  are  more  IFRA  than  H  .  Since  the  rank  tests  that  are  optimal 
for  the  Weibull  distribution 

*  b 

H(t)  -  1  -  exp(-at  ),  a,b  >0;  t  £  0 

coincides  with  the  optimal  tests  for  the  exponential  distribution  (Lehmann  (1953)), 
it  follows  that  the  tests  in  Sections  3  and  4  are  also  minimax  for  the  class  of  F 

A  A  U  1  A 

that  are  more  IFRA  than  H  .  Note  that  -log[l-H(t)  ]/t  ■  at  ,  thus  H  is 
IFRA  for  b  >  1  and  DFRA  for  0  <  b  <  1  . 
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7.  COMBINATION  OF  INDEPENDENT  TEST  STATISTICS.  GOODNESS  OF  FIT 

Let  T^,  . . .  ,  be  independent  and  such  that  has  the  continuous  distri¬ 

bution  function  G^*)  .  Let  F^(*).  •••»  F^*)  t>e  n  given  continuous  distribution 

functions,  then  the  problem  is  to  test  H  :  G.(*)  ■  F. (•).  i  *  1 . n,  against 

oil 

^1:  i  *  1,  n,  with  strict  inequality  at  least  once.  In  the 

goodness  of  fit  problem,  it  is  usually  assumed  that  F^  *  . . .  =  F^  ■  F  and 
G^  ■  ...  *  *  G  .  The  goodness  of  fit  problem  considered  is  not  the  one  in  which 

the  "shape"  of  the  distribution  is  of  interest,  but  where  one  is  interested  in  test¬ 
ing  whether  a  new  procedure  leads  to  stochastically  larger  variables  T^,  ..., 
than  a  standard  procedure  in  which  the  variables  have  known  distributions 
F^,  ...»  F^  .  For  references  to  the  "combination  of  independent  test  statistics" 
problem,  see  van  Zwet  and  Oosterhoff  (1967). 

A  test  <J>  is  called  monotone  if 


(7.1) 


tj^  >_  t^,  i  ■  1,  ....  n,  implies 
$(t^,  •••»  t^)  f.  $(t^»  •••»  in) 


and  it  is  called  strongly  distribution- free  (SDF)  if  its  distribution  depends  on 
^Fi’Gi^  only  through  G^^  *(•),  where  G^F^  *(t)  ■  P(Fi(Ti)  ±  t^)  •  Let 
6 (<^ ; F; A)  denote  the  power  of  the  test  $  for  the  IFRA  scale  alternative  in  which 
Gi(0  -  Fi(*/A1),  and  F^  )  is  an  IFRA  distribution.  Let  K±(x)  -  1  -  exp(-x), 
x  _>  0,  and  let  A  <_  1  (A  _>  1)  mean  A^  <_  l  (A^  1) ,  i  *  1,  ...»  n,  then 


Theorem  7.1 

If  <J>  is  a  monotone  strongly  distribution- free  test ,  then 

(i)  for  each  A  <_  1,  3 (d> ; F; A)  _<  B($;K;A),  and 

(ii)  for  each  A  _>  l,  e(4»;K;A)  j<  B(4>;F;A) 
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Proof : 

Set  «■  -logfl-F^(T^) ]  .  Then  the  SDF  property  implies  that  0($;F;A)  equals 
the  power  for  testing  H^:  "T^  has  the  exponential  distribution  K"  vs  Hj:  "T^ 
has  the  distribution  ( * )  -  F^F^  ^[K(*)]/A}"  .  Using  the  transformation 
T^'  *  -  A  log[l-F(T^/A)  ]  ,  the  IFRA  property,  and  the  arguments  of  Section  2,  one 
finds  that  G*(*)  <_  K(«/Ai)  whenever  A^  >_  1,  i  -  1,  n  .  Since  monotone  SDF 

tests  have  monotone  power  (Chapman  (1958)),  (ii),  follows.  The  proof  for  A^  _<  1, 
i  •  1,  n  is  similar. 

Fisher's  test  <J>  which  rejects  for  large  values  of  the  statistic 
r 

n 

(7.2)  -I  log[l-F.(T.)] 

i-1  1 

is  strongly  distribution-free  and  monotone.  It  is  also  uniformly  most  powerful  for 

1/0 

the  Lehmann  alternatives  G^(t)  ■  1  -  [l-F^(t)]  ,  9  >  1  (which  includes  the  ex¬ 

ponential  and  Weibull  scale  alternatives).  It  will  next  be  shown  that  it  is  minimax 

for  the  scale  model  when  A,  "...  -  A 

1  n 

Theorem  7.2 

(i)  Fisher '8  test  $  is  minimax  for  the  IFRA  scale  alternative,  i.e.,  if 

r 

each  F^  ranges  over  the  class  of  all  continuous  IFRA  distributions ,  then  whenever 

A,  *  . . .  ■  A  >  1  , 

1  n 

(7.3)  inf  B($;F;A)  <  inf  6($p;F;A) 
for  all  level  a  tests  <p  . 

(ii)  If  <j>'  is  another  level  a  SDF  test  satisfying  (7.3),  then  <fr'  ■  $p  a.e. 
Proof : 

Theorem  7.1  implies  that 

inf  8(*f;F;A)  -  0(<frp;K;A)  . 


(7. A) 
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Using  the  Neyman-Pearson  Lemma,  one  finds  that 
<7-5>  8($;K;A)  £  8(*f;K;A) 

for  all  $  such  that  6 (<#> ; K;  1)  £  a  .  Since 
<7-6)  inf  8($;F;A)  £  B($;K;A), 

the  result  (7.3)  follows.  The  uniqueness  part  follows  from  the  uniqueness  part  of 
the  Neyman-Pearson  Lemma  and  the  structure  (d)  property  of  SDF  statistics  derived 
by  Birnbaum  and  Rubin  (1954)  and  Bell  (1960).  The  structure  (d)  result  implies 
that  all  SDF  statistics  are  functions  of  T^,  ...,  T^  only  through  F^T^) . 

F(T)  • 

n  n 

Remarks 

(i)  Note  that  the  minimax  property  is  uniform  in  A  «=  A^  ■  . . .  *  A^  and 
that  it  holds  in  the  class  of  all  tests  satisfying  8(4>;K;1)  £  a  (not  just  the 
class  of  SDF  level  a  tests). 

(ii)  If  the  A's  are  not  all  equal,  then  the  minimax  statistic  is 
n 

-  I  [(A  -1)/A.)  log  [1-F . (T  ) ]  . 
i«=l  1  1  11 

A  sequential  test  is  monotone  if  t^  £  i  “  1»  2,  ...»  implies 

♦  (tj,  1 2>  ...)  £  4»(t|,  t^,  ...)  .  SDF  tests  are  defined  as  before.  Theorem  7.1 

holds  also  for  monotone  SDF  sequential  tests.  One  such  test  is  the  test  $  which 

1 

at  the  nth  stage  rejects  Hq  if  >  log  b,  accepts  if  T^  <  log  a,  and 

continuous  sampling  of  log  a  £  T  <  log  b ,  where 

r 

T  -  -[(A.-D/AJ  l  [1-F.  (T.)  ] 

*  i-1  11 

-n  log  A^,  A^  >  1,  and  a  <  1  <  b  . 


(7.7) 
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it 


This  is  the  sequential  probability  ratio  test  for  the  Lehmann  alternatives,  more¬ 
over,  it  coincides  with  the  sequential  probability  ratio  test  for  exponential  scale 
alternatives.  Let  E [ - | F , A ]  denote  the  expected  value  when  F^(*)»  F2(‘)i  ...  are 
the  Hq  distributions  of  T^,  Tj,  ...  and  F^(*/A),  FjO/A),  ...  are  the 
distributions.  If  N  denotes  the  number  of  stages  before  termination,  then  the 
optimality  property  of  sequential  probability  ratio  tests  and  the  above  implies  that 

Theorem  7.3 

If  a  and  b  are  ouch  that  E[$  | K ,  1  ]  ■  a  and  E[$  ( K , A .  ]  ■  8,  then  <J> 

A1  1  1 

minimizes  E [ N | K ,  1  ]  and  E[n|k,A^]  in  the  class  of  all  tests  $  satisfying 
E [ 4>  I K,  1  ]  <_  a,  E[4>  |  K,  A1  ]  _>  8,  E [N | K,  1  ]  <  <*>,  and  E [N  J K, Ax ]  <  »  .  Moreover ,  when 

A  *“ 

F^(-),  F2(0,  •••  are  continuous  IFRA  distributions ,  then  E [ 4>^  |f,A]  <_  a  for  each 

A  .  l 

A  <_  1,  and  E[<|>  |f,A]  _>  8  for  each  A  _>  A.  . 

A,  i. 
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